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Abstract. We present a new, short, self-contained proof of localization 
properties of multi-dimensional continuum random Schrodinger opera- 
tors in the fluctuation boundary regime. Our method is based on the 
recent extension of the fractional moment method to continuum models 
in 0, but does not require the random potential to satisfy a cover- 
ing condition. Applications to random surface potentials and potentials 
with random displacements are included. 

1. Introduction 

1.1. Motivation. We are concerned here with proving localization proper- 
ties of multi-dimensional continuum random Schrodinger operators in the 
fluctuation boundary regime. 

Such results were first found via the method of multiscale analysis, which 
had been developed in the 80s to handle lattice models and was later ex- 
tended to the continuum (for a rather complete history and list of references 
on multiscale analysis see [3U| and for some of the more recent developments 

Later, the fractional moment method was developed [3] as an alternative 
approach to the same problem, also initially for lattice models. It leads to a 
stronger form of dynamical localization than multiscale analysis (see ^ |3] ) 
and has provided much shorter and more transparent proofs in the lattice 
case, for example fT3] , 

It was recently shown in ^2 that all the main features of the fractional mo- 
ment approach also apply to continuum random Schrodinger operators. This 
extension required substantial new input from operator theory and harmonic 
analysis. The paper [2| provides a framework of necessary and sufficient cri- 
teria for localization in terms of fractional moment bounds, which can be 
verified for a rather broad range of regimes. 

One of our goals here is to complement the general framework from |2j 
by focusing exclusively on presenting a short and self-contained proof of 
localization properties via fractional moments for one specific regime, where 
the technical effort remains minimal. 

For this we pick a fairly general setting we label the fluctuation boundary 
regime. This is described by a random Schrodinger operator of Anderson- 
type in L 2 (M. d ), where our approach allows for quite arbitrary background 
potentials and geometries of the random impurities, provided the ground 
state energy is induced by rare events (fluctuations) and therefore sensitive 



Date: February 7, 2008. 



1 



2 A. BOUTET DE MONVEL, S. NABOKO, P. STOLLMANN, AND G. STOLZ 



to changes in the random parameters. The goal is to prove localization in 
the vicinity of the bottom of the spectrum. Of course, various versions of 
the fluctuation boundary regime have been studied in many works and we 
borrowed the term from |27| . 

Another motivation for our work is that we want to extend the fractional 
moment method to situations where the random potential does not satisfy a 
covering condition, i.e. where the individual impurity potentials have small 
supports which do not cover all of M. d . This condition, which was required 
for the technical approach to the continuum found in [2j, is not natural 
in the fluctuation boundary regime and should not be needed there as has 
already been verified via multiscale analysis. Particularly interesting exam- 
ples are random surface potentials which act in a small portion of space only. 
Nevertheless, they lead to a fluctuation boundary by creating new "surface 
spectrum" below the "bulk spectrum" . 

In our main result, Theorem ^ below, the fluctuation boundary regime 
will be described in form of an abstract condition. For random surface 
potentials, which will be discussed as an application, this condition follows 
in an appropriate setting from a result proven in |24j in order to derive 
Lifshits tails. Another application concerns models with additional random 
displacements as were originally studied in |1(J| . 

Let us confess that we require absolutely continuous distribution of ran- 
dom couplings. While it might be possible to relax this to Holder continuous 
distribution (as has been done in the lattice case, e.g. [I]), the fractional 
moment method is so far less flexible in that respect than the multiscale 
technique. In particular, see the variant of multiscale analysis adapted to 
Bernoulli-Anderson models recently developed in jS] and an application of 
similar ideas to Poisson models announced in jllj . 

1.2. Results. Let us now describe our results in more detail after intro- 
ducing some notation: On W 1 we often consider the supremum norm \x\ := 
maxj = i ... d\ x i\ an d write 



for the d-dimensional cube with sidelength r centered at x. For an open set 
G C W 1 we denote the restriction of the Schrodinger operator H to L 2 (G) 
with Dirichlet boundary conditions by H G . In our results we assume d < 3 
and rely upon the following assumptions, which guarantee self-adjointness 
and lower semi-boundedness of all the Schrodinger operators appearing in 
this paper: 

(Al) The background potential Vq G Lf oc unii (R d ) is real-valued, H := 



(A2) The set X C where the random impurities are located, is uniformly 
discrete, i.e., inf{|a — (3\ : a ^ G 1} =: rj > 0. 

(A3) The random couplings ij a , a £ I, are independent random variables 
supported in [0, r/ max ] for some rj ma _ x > and with absolutely con- 
tinuous distribution of bounded density p a with a uniform bound 
sup Q HpqIIoo =: M p < oo. 
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The single site potentials U a , «£l satisfy 

CUXA ru (a) <U a < CuXA Ru ( a ) 

for all a with cu, Cu,rjj, Ru > independent of a. 

Vw{x) = y~]r] a (u})U a (x) 

and 

H := H(u) := H + V u in L 2 (R d ). 

The most important condition expresses the fact that the ground state en- 
ergy comes from those realizations of the potential that vanish on large sets: 

(A4) Denote E := inferno) < mf a(H(uj)) and let 

Hp ■= Ho + r/ max U a , 

the subscript F standing for full coupling. 

Assume that Eq is a fluctuation boundary in the sense that 

(i) Ep : = inf a (Hp) > E , and 

(ii) There is m G (0, 2) and L* such that for m d := 42 • d, all L > L* 
and 

P(ct(# Al(:c V) n [£ ,£o + £~ m ] / 0) < £~ md - 

By Xx we denote the characteristic function of the unit cube centered 
at x. In the following it is understood that Xx(H G — E — ie)~ 1 x y = if 
Ai(x) fl G or Ai(y) n G have measure zero. 

Our main result is 

Theorem 1. Let d < 3 and assume (A1)-(A4). Then there exist 5 > 0, 
0<s<l,/i>0 and C < oo such that for I := [Eq,Eq + 5], a// open sets 
G cR d and x, y G M d ; 

sup Edlx^^-^-ier^n^Ce-^-"!. (1) 

Be/,e>0 

Exponential decay of fractional moments of the resolvent as described by 
([T]) implies spectral and dynamical localization in the following sense: 

Theorem 2. Let d < 3, assume (Al)-(A4) and let L be given as in Theo- 
rem Then: 

(a) For all open sets G C M. d the spectrum of H in I is almost surely pure 
point with exponentially decaying eigenf unctions. 

(b) There are [i > and C < oo such that for all x,y G M. d and open 
G c R d , 

E(sup|| X ^(# G )P/(tf G )xJ) < Ce-^. (2) 

where the supremum is taken over all Borel measurable functions g which 
satisfy \g\ < 1 pointwise and Pi(H G ) is the spectral projection for H G onto 
I. 
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Dynamical localization should be considered as the special case g(X) = e ltx 
in (b), with the supremum taken over t£l. 

The proof of Theorem ^ is given in Section |2J This will be done by a self- 
contained presentation of a new version of the continuum fractional moment 
method. While we use many of the same ideas as |2j, due to the lack of a 
covering condition we can not rely any more on the concept of "averaging 
over local environments", heavily exploited in It is interesting to note 
that, in some sense, we instead use a global averaging procedure. Techni- 
cally, this actually leads to some simplifications compared to the method 
in |2j, as repeated commutator arguments can be replaced by simpler iter- 
ated resolvent identities. We also mention that exponential decay in Jj) will 
follow from smallness of the fractional moments at a suitable initial length 
scale (the localization length) via an abstract contraction property. 

As technical tools we need Combes-Thomas bounds (in operator norm as 
well as in Hilbert-Schmidt norm) and a weak-L 1 -type bound for the bound- 
ary values of resolvents of maximally dissipative operators, which is based 
on results from |26j and was also central to the argument in [2j. We collect 
these tools in an appendix. 

That Theorem |2] follows from Theorem ^ was essentially shown in [5] , 
Section 2. In Section |3] below we will briefly discuss the changes which arise 
due to our somewhat different set-up. In particular, the argument in [5] for 
proving (J2J) uses the covering condition 



in one occasion. But this is easily circumvented. 

In Sections 0] and |S] we apply our main result to concrete models by veri- 
fying assumption (A4) for these models. In Section 0] we consider Anderson- 
type random potentials supported in the vicinity of a lower-dimensional 
surface. The "usual" fully stationary Anderson model is considered in Sec- 
tion [5J The fact that we don't have to assume a covering condition leads 
to high flexibility in the geometry of the random scatterers. We could use 
this to go for far reaching generalizations of Anderson models. Instead, we 
restrict ourselves to the treatment of additional random displacements as 
was done in 

1.3. Remarks. We could have extended Theorem \I\m at least two differ- 
ent ways, but refrained from doing so to keep the proofs as transparent as 
possible: 

(i) The restriction to d < 3 is not necessary. We use it because in this 
case the abstract fractional moment bound in Corollary El is more 
directly applicable to our proof of Theorem ^ than in higher dimension 
(which technically can be traced back to the fact that Xx(— A + 1) _1 is 
a Hilbert-Schmidt operator only for d <3). In higher dimension more 
iterations of resolvent identities would be needed to yield the Hilbert- 
Schmidt multipliers required by Corollary El leading to more involved 
summations in the arguments of Sectional 

(ii) Instead of bounded U a we can work with relatively A-bounded U a , i.e. 
allow for suitable L p -type singularities in the single site potentials. In 




(3) 
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the course of our proofs they could be "absorbed" into resolvents using 
standard arguments from relative perturbation theory. 

In principle, our arguments can also be used to prove localization at fluc- 
tuation type band edges more general than the bottom of the spectrum 
without using a covering condition as in [2|. But this would require to be 
much more specific with settings and assumptions and, in particular, with 
the geometry of the impurity set. Inconvenience would also arise from having 
to work with boundary conditions other than Dirichlet. 

We mention that the applications in Section 0] improve the results on 
continuum random surface potentials of EI] , obtained through the use of 
multiscale analysis: 

(i) The exponentially decaying correlations of the time evolution, shown 
as a special case of Theorem I2Jb), are stronger than the dynamical 
bounds which follow from multiscale analysis. 

(ii) Due to the use of the recent result of Q1Z on Lifshitz tails for surface 
potentials, we do not need a condition on the smallness of the distribu- 
tion of the r\ a near the minimum of their support as in [Jj, a progress 
that had been achieved in |24j . 

(iii) We can allow for more flexibility concerning the geometry of the scat- 
terers. 

Of course, due to using fractional moments we cannot include single site 
measures as singular as the ones considered in but instead have to 
assume absolute continuity of the rj a . 

2. Localization near fluctuation boundaries 

This section is entirely devoted to the proof of Theorem 1. For a conve- 
nient normalization write 

£ a (u) := ry max - 7y a (w) 

for LO = (u a ) ae i = (Va(u))ael 6fi:= [0,7/ max ] X , 

and denote the product measure <8> a ez dr] a p a (r] a ) on Q by P. We write 
W(x) := W u {x) := J2U")U a (x). 

Note that W u > and that 

H = H{u) = H F - W u . 
Fixing an open set G C M. d we write 

R G = Rf = (H G -z)- 1 , 

Rp = Rp z = (Hp — z) 1 

whenever z = E+ie. Since Hp 1 > Hp due to our choice of Dirichlet boundary 
conditions, and Ep = inf a(Hp) we know that (— oo, Ep) C p(Hp). 
The resolvent equation yields 

rG = r g + rG wr g + rG wr g wr g^ ( 4 ) 



IS 

Ihs 
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an identity that will be used over and again. The other workhorse result is 
the following averaging estimate, that follows from Corollary 1171 in the ap- 
pendix below, taking into account the uniform boundedness of the densities 

Pa- 

Lemma 3. For all s £ [0, 1) there is c(s) such that 

dVaPa(Va) J d W (^) \\M X U X J 2 {H G - E - ie)- 1 U 1 J 2 M 2 \ 

<c( S )||M 1 ||^||M 2 ||^ s . 
As a warm-up, we prove boundedness of fractional moments: 
Lemma 4. Let E\ < E F , I = [E , Ex] and s G [0, 1). Then 

sup{E \\XxRE+ieXy\\ S I E £ I, e > 0, x, y £ R d , G C R d open} < oo . (5) 

Proof. We use (@J) above and write, suppressing the superscript G and the 
subscript z = E + ie mostly: 

XxRXy = XxRfXv + XxRfWRfXv + XxRfWRWRfXv 

The first two terms on the r.h.s. of this equation obey an exponential bound 
due to the Combes-Thomas estimate, see subsection IA. II below: 

\\XxRFXy\\<ce-^ x -y\ 

and 

\\XxRFWR F Xy\\ < r] max Y,\\XxRFU^ 2 \\ • \\UV 2 R FXV \\ 

ael 

with hq and p\ = po/2 depending on E\ only. In the last estimate we have 
used that X is uniformly discrete. 

For the third term, expand W = Yla^U a and use the boundedness of 
the £ a and the fact that 

(5>") S ^5>n 

to estimate 

WXxRfWHWRfXvW < c \\XxRFU a RU p R F Xy\\ s . 

a, pel 

We now fix a, /3 6 X and use the workhorse Lemma |3] to conclude 



dn a p a (n a ) J dn l3 p f3 (r]p)\\xxRFU a RUf 3 R F Xy\\ s 

I^C^WXxRFU^W^WU^RpXyWks 
< c(s) ■ e — */*ola= — of] . e -sna\y-i3\ 

by the HS-norm Combes-Thomas bound from Proposition ^] and since 
dist(x, supp U a ) > \x — a | — Ru where Rjj majorizes the size of the sup- 
port of U a according to assumption (A3). 

Note that here and in the following we use the convention that c, c(s), 
etc. denote constants that only depend on non-crucial quantities and may 



LOCALIZATION NEAR FLUCTUATION BOUNDARIES 



7 



change from line to line. In particular, the constants are independent of 
e > and the random background. 

Now, we can sum up the last terms and get the assertion. □ 

Remarks, (i) In this proof it is still quite easy to see how to extend to 
arbitrary dimension through iterations of the resolvent identity. It will be 
harder to keep track of this later. 

(ii) Note that due to the a, /3-summations, averaging over the rj a is re- 
quired for all a, i.e., is global. In [2], due to the covering condition, an 
argument is provided that only requires averaging over local environments 
of x and y and proves Lemma 0] for arbitrary finite intervals / = [Eq,Ei], 
i.e. without requiring E\ < Ep. 

(hi) The above proof shows that © also holds in HS-norm, but this will 
not be used below. 

We will now start an iterative procedure that will show exponential decay 
of E(||xz-RXyl| s ) m \ x ~ y\ f° r energies sufficiently close to Eq. Clearly, it 
suffices to consider x,y £ 1j d . In view of the preceding lemma the following 
quantity is finite: 

T x>y := sup{E \\XxRE+ieXy\\ s \ E £ I, e > and G C JR d open} . 

Moreover, we should actually keep in mind the dependence on the interval 
/ = [-Bo, E\]. In fact, E\ will later be chosen small enough. 

In order to use that Eq appears rarely as an eigenvalue for boxes of side 
length L we exploit the resolvent identity and what is sometimes called the 
Simon-Lieb inequality in a way visualized in Figure l! 




Figure 1. The geometry of the induction step 
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Consider 

B L := A L (x)nG, 

dB L = (A L (x) \ A L _ 2 (x)) HG, and x L ■= XdB L - 
Furthermore, with Ru as in assumption (A3), define 
C L :=G\ A 2 r u+L (x), 

dC L = (A 2Ru+L+2 (x) \ A 2 r u+ l(x) ) n G, and \\ '■= XdC L - 

The geometry is chosen in such a way that R Bl and R L are stochastically 
independent. For R Bl we can use the fluctuation boundary assumption to 
get small fractional moments and the right size of L will be adjusted. But 
all that later... 

Thus, by the Simon-Lieb inequality (e.g. |H0] , Sect. 2.5) 

\\x x R G Xy\\ < c\\ Xx R Bl xI\\ ■ \\xIR G xt\\ ■ \\xtR CL Xy\\ (SLi) 

where C only depends on sup|^ Q | ag j} ||V ||oo and the interval /. 

The basic idea for proving exponential decay of T x>y is to establish a 
recurrence inequality for energies sufficiently close to Eq. This recurrence 
inequality is described in Proposition below and allows to apply a discrete 
Gronwall-type argument found in Lemma below. To this end we exploit 
smallness of fractional moments of the first factor on the r.h.s. of (SLI) for 
energies close to Eq and sufficiently large, but fixed, L: This will follow from 
(A4)(ii) as is presented in the following Lemma El Fractional moments of 
the second factor are bounded due to Lemma |U (up to a polynomial factor 
in L). Finally, we use the third factor to start an iteration (with x replaced 
by sites x' covering the layer Al + r u+2 \Al + r v ). By construction, the first 
and third factor on the r.h.s. of (SLI) are probabilistically independent. 
Unfortunately, the second factor introduces a correlation which prevents us 
from simply factoring the expectation. We will rely on a version of the re- 
sampling procedure developed in [2] to solve this problem. Moreover, we will 
not use Lemma HJ but apply Lemma El directly to bound certain conditional 
expectations. This will result in Proposition H3 below. 

Lemma 5. For m as in (A4) and s £ (0, |) there is L* = L*(m,s) such 
that for all L > L* , open B C A L (x), E € I := [E ,E + L~^ m ], e > and 
u, v € Z d with \u — v\ > -j we have 

n\\Xu{H B - E - ieT^W*) < L-l™*, 

where m& = 42 • d. 

Proof. Divide £1 into the good and bad sets 

^good := | dist(a{H B ), E ) > L~ m }, fl had = O \ fl good . 

Since H B > H^ L ^ by our choice of Dirichlet boundary conditions, (A4) 
implies that 

P(tt bad ) < L~ md . 

We split the expectation into contributions from the good and bad sets: By 
the improved Combes-Thomas bound Subsection IA.1I we get, for lj 6 r2g OOC i, 
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Eel: 

IIXu-flB+ieXull S CJv2 e 

This gives a uniform bound of the same type for the expectation over f2 g0 od- 
For the bad set Holder with t € (s, 1) gives 

H\\XuR% +i£ Xv\\ s Xn b J < (E(|| XuJ R| +i£ x,||*)) f P(^bad)^ f 

<c(t)-tL^-^. 

Now we choose t = |(s + 1) so that (1 — |) > ^ if s < ^. Putting things 
together, we get 

n\\xuR% +i£ xv\n < c(i)(LH- a H-i-^ + i(HH). 

If L is large enough we can use < 1 and \u — v\ > h to see that the r.h.s. 
is bounded as asserted. □ 

The exponential decay of the T xy will follow from the following result, 
whose proof will take most of the present section. 

Proposition 6. There exist L* , k > 0, c > and C > 0, all depending 
on s, m, Ru,ru, M p , Eq, Ep, ry max , such that for L > L* and I = [Eq, Eq + 
L _ 2 m ] the above defined t XjV satisfy: 

r x , y <L- 2d ~' i e-<\ x - x '\+\y-y'\)/ L tjj + Ce- c \ x ~ y \' L . (6) 

x',y'eZ d 

Proof of Proposition We now restrict to the energy interval / = [Eq, Eq + 
L~2 m ] assuming L is large enough to guarantee that I C [Eq,Ef). Using 
(SLI) above and denoting 



T X ,L 


= XxR Bl x l 


S X ,L 


= xlR G xt 


Qx,L 


= x1R Cl x v 



we get that 

n\\XxR G Xy\\ s ) < CE(\\T X , L \\ S \\S X , L \\ S \\Q X , L \\ S ). 
Note that ||T XjL || s and \\Q 

X) l\\ s are stochastically independent. Unfortu- 
nately, they are correlated via \\Sx ,l \\ 8 • 

Fix s € (0, g) to estimate E(||T Xi £|| s ). Using the preceding Lemma, we 
get that 

n\\TxA s ) < E n\\xxR BL x z \\ s ) 

< CL d - 1 -L-^ md , 
for L large enough. We get that 

E(||T X , L || s ) < L d -^. 
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We can now expand S x ,l to split off a uniformly bounded (in lo) term: 

s x , L = xlRfxi + xIR g f wr g fX + + x - l R$wr g wr g fX + l ■ 



Since I C [Eq,Ef) we have that HS^lII 5 is uniformly bounded. Thus, we 

get 

m\X,R%eXy\\ s ) < C{E(\\T XjL \\ s ■ \\Q X , L \\ S ) + S 2 ) 

= C{E(\\T XtL \\ s )-E(\\Q x , L \\ s ) + Z 2 ) 

as \\T x< l \\ s and are independent. Here 

S 2 -M(\\T X>L \\ S \\S 2 , L \\ S Q X>L \\ S ). 

Expanding x£ we S e t> f° r some c > 0, that 

x'£dC L 

whence 

T,,„ < Z, J -i"" V e-*- I 'l/ 1 T J ,.„ + C sup E 2 . (7) 

To estimate £ 2 we begin by expanding 
T X , L = XxR Bl Xl 

= XxR F L Xl + XxR B F L WR B F L xl + X x R B f L WR Bl WR B f l Xl - (8) 

Since I has positive distance from a(H F ), we have the Combes-Thomas 
bound Ce _W)L / 2 for the norm of the first two terms on the r.h.s. of (JSJ), see 
Appendix IA.1I Here C < oo and > are uniform in the randomness, 
E G I, £ > and x € Expanding the third term and using boundedness 
of the £'s yields 

ii^,l|| s <c(V^-*+ iit A7 i 

/3,~,£inA L+Ru (x) 

setting 1>, 7 = XxR^R^U.R^xl, and only summing over those /3, 7 
for which the corresponding [/-terms touch B F . 
A similar argument applied to Q x ^l leads to 

\\Qx,l\\ s < c(e-">"<.\*-v\-%) + Yl ll^',7'll s ), 

where we have chosen Q/j'y = XLRF L U-y'R Ch UpiR F L Xy 
Finally, expand 

S 2 ,l = x~ l RfWR g WR$x + l = Y, 
where S a ^ = xiR G £ a U a R G £, a >U a >R ( F 'xi- 
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Combining all this we get that 



^ - ■ us Q -no-s-{\x-y\-^) 



a,a' 
L 



+ e-w-^ Y n\\s a ,A 8 -\\Qp 

+ Y EdlT^IMIS^ine 



Ho-s-{\x-y\-^) 



a,a' ,/3,7 

+ £ e(||t A7 |M|s q ^|I s - IIQ^7'II S ))- 

a,a',/3,/3',7, 7 ' ' 

The most complicated of these terms is the last one; it will be obvious how 
to estimate the first three once we have established a bound for the last one 
according to the assertion of Proposition E3 Thus we have to estimate 

S3 := Y A a,a',/3,P',j,Y, 
a,a' ,/?,/?', 7, 7' 

where 

■A a , a ',p,p',j,j' = E (||T/3 i7 || s • ||£ a)Q /|| s • ||Q/3', 7 '|| S ) • 
If it weren't for the »S QjQ ,'-terms, the Tp n and Qp^y would be independent, 
leading to an estimate like in (JJJ) above. We will reinforce a certain kind of 
independence through re-sampling. For fixed 

J := {a, a', 7,7'} 

we introduce new independent random variables £j, j € J, independent 
of the £c> C e X, and with the same distribution as the We denote 
the corresponding space by Q, the corresponding probability by P and the 
expectation with respect to P by E. Consider 

H(u,Q) = H{u) + Y&{o>) ~ U")) u i 

" v ' 

w 

and note that H doesn't depend on the £j, j E J . The resolvent identity 
for Rf = (H G — z)~ x gives 

Rf = + RfWRf. 
We insert this for T^ j7 and Qp* n > and get 

T Pn = XxRpUpR^U^xl+XxR^UpR^WR^UyRpxl 
s * ' ' ' 

T /3,7 T /3,7 

and, similarly, 

Q(3',y = Qp'rf + Q/3',7'- 

Now we can estimate 

Aa,a.',p,p'aa' <EE (||T^ )7 || S + ||T/3, T ||*)||i?a,a'[| s (||Q^7'|| a + [iQ^.yll*) ■ 

(9) 
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This gives a sum of four terms we have to control. Let's start with the 
easiest one 



Denote 



A i,a',l3,f3', 7,7' : = E E [ll^/3, 7 iri|5 , Q ! ,a'iril^/8',7'l| S ] • 

E(X\a,a') = J d£aPa(U J d£cfPa>{Zo)X{Z) 

for a random variable on Q x ft, so that E(X|a,a') is nothing but the 
conditional expectation with respect to the cr-field generated by the family 
| (3 G X \ {a, a'}). According to the usual rules for conditional expecta- 
tions: 



= EE 



E(||r^lS Q , a nilQ^y||>,a') 
\\TpJ s \\Q p , jY \\ s E(\\S a , a ,\\ s \a,a'} 



(10) 



since the T and Q are independent of £ a , £ a /. Using the workhorse Lemma El 
and the Combes-Thomas estimate Proposition ^] we get 

E(\\Sa,a>W>> a ') < c(s)||xZ^i||Hsll^^FXillHS 

< c ( s )L 2s ( d - 1 ) e - Atls(| t- |a - x|l + | t- |a '- :E|l) , 

where the extra L 2s ( rf_1 ) term comes from covering 8Bl and 8Cl- We have 



EE 



f/j, 7 ||lQ/3' l7 '|| s =E[||T A7 || s ||Q^y|| s ] = E[||7> i7 || s ]E[||Q /3 ,. 



since the £'s have the same distribution as the £'s and the T's and Q's are 
independent. Inserting into IjlflJI gives 

< Q ',A/3',7,y < c^L^-^e-^^lt-l^ll+lf-l^-IDEfllr^inEfllQ^^H 8 ; 
We will now treat the latter two terms separately: 



Step 1. Denote by Zi^f) = {y' G Z d | Xy' ' / 0} those lattice points 
whose 1-cubes support t/y . By Combes-Thomas once again: 

WQytfW = \\xiRF L U Y R CL U p/ R^ X y\\ s 

<C Yl \\xx>R CL Xy4 s z~^ x ~ A ~ L) e~ t * lAy ~ y '\- 

x>£Z(y)yi£Z{f3>) 

By assumption on the size of the support of Uy we see that #Z("f') is 
uniformly bounded. This and uniform discreteness of I gives 

^EHQ^yH 5 < C e~^ x - x ' l -^e~^ y - y ' l T x ^ yl . 
P'rf x', y 'ez d nc L 
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Step 2. For the T^ i7 -term we have 



IT, 



A7I 



\\ Xx R% l U l3 R BL U 1 R§Pxl\\ a 
^ C E E Wx,RF L Xu\\ s \\xuR BL Xv\\ s \\xvRf L xl\\ s - 

ueZ(/3)r\B L v&z(j)r\B L 
If \u — v\ > jL, Lemma gives 

n\\XuR BL Xv\\ s ) <C-L-* m *. 

If, on the other hand, \u — v\ < |L then dist(u, <9-Bl) > or |a; — u\ > |L, 
so that the uniform bound of Lemma 0] for E(||xui? Xu|| s ) together with 



the Combes-Thomas bound for \\x v Rf L Xl ll s > resp. HXx-R^Xull 8 gives, for 
L large enough, 

^{\\x,Rf L Xu\\ s \\XuR Bl Xv\\ s \\XvRf L xI\\ s ) < Ce~^ sL 



< L~2 



Combined we get that, again for L sufficiently large, 

Y J n\\T Pn \\ s )<cL id - 1 2 m \ 



*,7 



where an extra factor L 2d arises through the number of terms considered. 
Joining Step 1, Step 2 and the bound 



^ e -Mis(i-Ha-z|i+H-K-z|i) < C ^ L 



2d 



we arrive at 



E 4 



a,a',/3, / 3',7,7 / 



< C(s)L 6d -^ md ^ 

x',y'eZ d nC L 



e -^is(\x-x'\-§) e -fns\y-y'\ ; 

x ,y j 



which is a contribution to S3 (and therefore S2) bounded by one of the type 
asserted in Proposition H3 

A look back at © shows that we still have to estimate three terms similar 



to ^a,a',/3,/3',7,7' °f wm ch the last one, 



A a,ce', /3,/3', 7,7' 



EE 



|T/3, 7 || S ||«S'a,a'ini < 3/3' 



is the most complicated one. Using Steps 1 and 2 above as well as the steps 
below it will be clear how to treat the two remaining terms. 



Step 3. We start taking the conditional expectation: 

E(||7> i7 |n|5 Q , a / W'WQprf \\ s \a, a', 7, 7') 



A i,a',/3,f3',-y,Y ~ E E 
< EE 



EfHT, 



3s 



I3 r/ W |a,a,7>7J 3 
• E(||S a , a , || 3s |a, a', 7, y)3 • EfllQ^y || 3 >, a', 7, Y)* 
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by Holder's inequality. Like above, the middle term can, up to C L 2s<yd ~ l \ 
be estimated by 

Recall that 

\\Qfy,T W 3s = \\x!Rf l UyR Cl Y & " ^) u ^ CLU P' R F L x y \\ 3s 

<C- Y WxlR^UyR^U^U^R^Xyf 3 , 
jeJ\{~t} 

where 7 can be excluded from the summation as t/ 7 doesn't touch Cl- 
Integration over £j and £y gives a uniform bound by the workhorse LemmalHl 

E(\\Qp',7'\\ s >,a',^') 

< Y mxlR^UyR^UjR^UpR^Xyf^a'nn') 

jej\{ 7 } 

<C(s)- Y \\xlR C F L U^ s -\\u}R^U p/ R c F L Xy \\li s 

jej\{-y} 

so that, as the sum has only three terms, 

E(||Q^ y || 3 >,a', 7 , 7 ')s <C7( S )- Y ||^ l ^IIhs • \\uf R c ^R^ X y\\ s us ■ 

j&J\{-y} 



Similarly, 

E(||f /9)7 || 3 '|a )Q ',7 ) y)i < C(s)- Y Wx,RF L Uf3R BL Uf\\ s m ■ R^xlWks 

j&J\{i'} 



Now and Eq are independent so that 

<«W 7 ,y < C(s)L 2s ^ ■ EE[E T ] • EE[S Q ] • e a>Ql . (11) 

Since the £j and the £j have the same distribution, we can omit the hats in 
R Cl and R Bl and replace EE by E in the bounds for EE[S T ] and EE[Sq] 
to be derived below. 

Step 4- We start with the Q-term. Combes-Thomas, Proposition 1131 gives 
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This will be used to deal with the term for j = 7' which appears in the sum 
over J~\ {7}; since ||AB||hs < \\A\\ II-^IIhs we get that: 

IIx^C^IIhS • \\U^R CL U pl R^ X y\\ S HS 

< CL< d -^%-^ s ^'- x \-^\U^R CL ul\\ s ■ \\Uj,R c F L x y \\ks 
< CL »(d-i) ^ e-^ x - x ' hlMlsly - y ' l \\xx'R CL Xy'\\ s - (12) 

For the terms j = a and j = a' in the sum we borrow from e a ^ a i above and 
use that 

e 3 f < C . e -c\x-x'\/L 



if j € {a, a'} and x' £ Z(j): 

4,JxlR^\\ s ■ \\ufR c ^R C F "xy\\ks 

< CL 8(d-l) e -n,\\>y'-x\-% e-^ X - X '\-^ y - y '\\\x x 'R CL Xy'\\ S - (13) 

y'ez(p') 

1 

Summing each of the three contributions from (|12|) and to ea jQ XQ over 
/S',7 7 (and extending the x'-sum in (fl3|) to all of Z d ) gives 

£ e J Q ,EE[E Q ]<C'L 2 ( d - 1 > e-rl— 'l-M lS |y-y'l Tl , y . (14) 

/3',7'ex x'^'eZ 1 * 

We now show that summation over a, /?' gives a small prefactor. By 
exponential decay: 

E e L'<^ M (15) 

a, a' 

.Step 5. We analyze 

\\E( Xx R^UfsR^u}\\ s HS ) < llx^^lll&s-ECIlC^^^/n. 
If \I3 - j\ < \L then either \x - (3\ > \L or dist(j,<9C L ) > \L. Since 

either EiWxxR^UpR^Uj \\ s ns ), e\ a , or \\U* B*P xlWhl 

is bounded by L~2 md ; see Step 2 above. If, on the other hand \f3 — j\ > \L 
we can use Lemma El above to estimate 

E(\\UjR BL Uj\\ s ) <C-L-^ md . 
Summing up these terms we get that 

]T e| a , EE[£ T ] < CL 3 *- 1 -*™* (16) 
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since (3, 7 run through at most cL d different points of X in B^. Putting the 
estimates from (|14 p . (jl6| ). (|15j) together we arrive at: 

a,a',P,P',%-y' x',y'£Z d 

which is the desired bound. To deal with the other terms appearing in 
A^a',/3,/3',7,7' we just combine the corresponding steps to control the T and 
Q-sums respectively. 

This concludes the proof of Proposition El □ 

For energies sufficiently close to Eq we will now complete the proof of 
exponential decay of T x>y , and thus of Theorem ^ by applying a discrete 
Gronwall-type argument to the recursion inequality established in Proposi- 
tion 

For \x > consider the weighted £°° -space 
i.e., for if) = (i> x>y ), 

U\\ x = sup ^ x -y\i 2 \^ y \. 

x,y£Z d 

Lemma 7. The operator A defined by 

x',y' 

is bounded as an operator in X as well as an operator in with 

\\A\\ X < C(d)fi~ 2d and \\A\\t»o < C{d)^ 2d . (17) 

Proof of Lemma ^ The norm of A in X is the same as the norm of the 
operator A in ^°°(Z 2d ) with kernel 

2, , , - e lAx-y\/2 -n{\x-x'\+\y-y'\) -iAx'-y'\/2 

Thus 



\A\x = ll^lk 00 = sup A xyx i y i 

x,y , , 
x',y' 

< csup [[ dx'dy'e^-y^e-^-^y-y'^e-^'-y'^ 2 



x,y J J 

= Csup J J ds dpe^ A l/ 2 e-^ s l + lf- A D-^- s l/ 2 , 

with the substitutions s = x — x',p = y' — x, A = y — x. 
Bound the latter exponent through 

fi(\s\ + \p-A\) + ^\ P -s\ 



(m - ~ )(H + \P~ A|) + |(|A -p\ + \p-s\ + \s\ 



>|(| s | + |p-A|) + ||A|. 

After cancellation the integral factorizes and gives ()17|) for ||j4||x after scal- 
ing. The bound for is found more directly. □ 
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This may be applied to the situation of Proposition as it shows that for 
L sufficiently large the operator A with kernel 

A , , _ T-2d-K -c(\x-x'\+\y-y'\)/L 
-^xyx y ^ c 

has norm less than one, both as an operator in X = £°°(Z 2rf ; e c \ x ~y\/ 2L ^ and 
an operator in £ OD (7 J 2d ). Fix this L and choose 5 = L~ m , I = [Eo, Eq + S] in 
Theorem n and the definition of t X)V . 

The recursion inequality © now takes the form 

r x , y < (Ar) Xt y + b x>y (18) 

with b x>y := Ce~ c \ x ~ v ^ L . The conclusion of the proof of Theorem ^ is now 
the content of 

Lemma 8. 

t = (r Xj y) G X. 

Proof of Lemma\$ With /i = define the diagonal operator 

V = diag^ 1 ^/ 2 ), 

which is an isometry from X to l 00 {l? d ). Let f = Vt and b = Vb G 
Let A = T>AT>~ 1 . Then (|18j) implies that componentwise 

f<Ar + b. (19) 

Since r = T>~ l f is bounded and A a bounded operator in l ac '(T? d ) : we have 
that f G Y := r°(Z 2d ; e^l^l/ 2 ) and i is a bounded operator in F with 
non-negative kernel. Thus we obtain from (|19|) that 

i"f < i™ +i f + ah 

holds with finite components. Summation yields 

JV 



f <A N+1 r + ^A"6 



n=0 

and thus 

N 



T <A N+1 T + Y,A n b 



n=0 

for all N. 

A: e°° -> £°° is a contraction and r G r°. Thus A^+H -» in £°° 
and componentwise. Also, A: X — » X is a contraction and b £ X. Thus 
X^=o ~~ * ~~ A) _1 6 G X and componentwise as JV — » oo. We conclude 

r < (I-A)~ l b G X. 

Lemma |S1 is proved. □ 
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3. On the proof of Theorem [21 

That the localization properties stated in Theorem [2] follow from the 
fractional moment bound for the resolvent established in Theorem ^ was 
demonstrated in (2j . Here we want to comment on two minor changes in the 
argument which are due to our somewhat different set-up. 

First we note that spectral and dynamical localization as established in 
parts (a) and (b) of Theorem El hold for restrictions of H to arbitrary open 
domains G, and, in particular, that the exponential decay established in 
equation (J2J holds with respect to the standard distance \x — y\ rather than 
the domain adapted distance distG*(x, y) used in 0. Given that the corre- 
sponding bound Q in Theorem ^ is true for arbitrary G and in standard 
distance, this follows with exactly the same proof as in Section 2 of (with 
one exception discussed below). That the authors of [2j chose to work with 
the domain adapted distance was in order to include more general regimes 
in which extended surface states might exist. This is not the case in the 
regime considered here. 

Second, let us provide a few details on how to eliminate the use of the 
covering condition |2J) from the proof of (J2J) provided in Section 2 of [21 . As 
done there one first considers bounded open Acl d and defines 

Y A (I;x,y):= sup \\ Xx f(H A ) X y \\- 

/6C C (/),|/|<1 

If E n and ip n are the eigenvalues and corresponding eigenfunctions of H A 
and / is as above, then f(H A ) = Yl n -.E n ei f( E n)(ipn, • )4>n readily implies 

Y A (I;x,y)< ^2 \\Xxipn\\ ■ \\Xy^n\\- 

n:E n £l 

At this point we modify the argument of [2j and write 

Xy^n = Xy(H F - E n )-\H A - E n )^ n 
= Xy {H£ - E n )~ l W^ n 
= Y J ^Xy(H A -E n )- 1 U a ij n . 

As all E n € / have a uniform distance from inf a (Hp) we get from Combes- 
Thomas Proposition El that 

llx^nll < Cj2\\Xy(H F - E^U^W ■ \\Ul^ n \\ 

a 

<cY,^ My ~ a{ \\ul /2 M. 

a 

Inserting above yields 

Y A (I;x,y) < Cj2^ o]y ~ a] Qi(I;x,a) 

a 

with Qi(I;x,a) = Yl n -E el llXx^nll ' \\U a ipn\\ defined as in [2J, where the 
bound E(Qi(J; x, a)) < C'e - ' 11 ^ - "' is established without any further refer- 
ences to the covering condition. Thus we conclude 

E(Y A (I;x,y)) <Ce-^ x -v\. (20) 
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The rest of the proof of Theorem |21 in particular the extension of (|2U[) to 
infinite volume and a supremum over arbitrary Borel functions, follows the 
argument in j2] without change. 

4. Localization for continuum random surface models 

Random surface models have attracted quite some interest with most of 
the work dealing with the discrete case El El El El El E] and some 
with the continuum case |2L)l 1241 as we do here. Our aim in this 
section is to show that under suitable conditions such surface models obey 
condition (A4) above. To achieve it, we combine recent results from j^H 
with a technique from 29 . 

As usual, the background is assumed to be partially periodic: 

(Bl) Fix 1 < di < d and write R d = R dl x R d2 , x = (x 1 ,x 2 ); assume that 
Vq € Lf QC unif is real- valued and periodic with respect to the first 
variable, i.e., 

Vq{x\ + m, X2) = Vq(xi, X2) for m 6 Z dl . 

Denote H := -A + V . 

In order to state our second requirement, let us recall some facts from Bloch 
theory. For more details, see [21]. For Vo,Hq as in (Bl) we get a direct 
integral decomposition 




where T dl = M rfl /(27rZ) dl is the d± -dimensional torus and 

h e = -A + V in L 2 (5i) 

with ^-periodic boundary conditions on the unit strip S\ = Ai(0) x W 12 . We 
now fix the assumption 

(B2) 

inf a(ho) < inf a ess (h ). 
It is well known that under (B2) we have that 

Eq := inf a(Ho) = inf <r(/io) 
and there is a positive eigensolution ipQ of the distributional equation 

H tp = E ip , 

see |241 123j and the references therein. Finally, our random perturbation is 
assumed to satisfy 

(B3) The set I C M. d , where the random impurities are located, is uniformly 
discrete, i.e., inf{|a — j3\ : a ^ G 1} =: rj > 0. Moreover I is dense 
near the surface M. dl x {0} in the sense that there exist R± , c± > 
such that for L large enough and x\ E Mr 1 : 

#[lfl(A L (n)xA fli (0))] >c L L d \ 

We will see that (B1)-(B3) ensure (A4) from Section 1. Of course, there 
might be other ways to verify (A4) for surface-like potentials so that Theo- 
rems 1 and 2 could, in principle, be used for other examples. 
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Theorem 9. Assume (B1)-(B3) and (A3). Then there exist 5 > 0, < s < 
1, \i > and C < oo such that for I := [Eq, Eq + 5], all open sets Gc l d 
and x, y G M. d 

sup n\\x,{H G -E-\e)- l Xy f)<Ce-^y\. (21) 

Eel,e>0 

In particular, the following consequences hold: 

(a) The spectrum of H G in I is almost surely pure point with exponentially 
decaying eigenf unctions. 

(b) There are n > and C < oo such that for all x,y € Z d , 

E( S up\\ Xx e- itHG Pi(H G )xy\\) < Ce-^. (22) 

teffi 

The rest of this section is devoted to deducing (A4) under the assumptions 
of the Theorem. Note that this will be accomplished once we have shown 
the following, where 

S L = S L ( Xl ) := A L ( Xl ) x M d2 

denotes the strip of side length L centered at x± € M. dl perpendicular to the 
"surface" R dl x {0}. 

Proposition 10. For all 7, £ > there exists L(j,£) such that for all odd 
integers L > L(7,£) and x\ € 

F{a{H s ^) n [Eq, E + IT*] + 0} < IT*. (23) 

In fact, (A4)(ii) then follows, since 

H k l{x) > H S L ( Xl ) and therefore E Q < 

\nia{H s ^) < \nia{H K ^). 

We will actually prove the analogue of Proposition El] with Dirichlet 
boundary conditions replaced by suitable Robin boundary conditions that 
are defined using the periodic ground state ipo introduced above. Assume, 
for later convenience, that 

/ |V>o(*)| 2 dx = l. 

J Si 

We consider on S^, L £ 2N — 1, Mezincescu boundary conditions, given as 
follows. Let 

X(x) ■= — ——-V n ^o(x), 
ip (x) 

where V n denotes the outer normal derivatives. The Mezincescu boundary 
condition can be thought of as the following requirement for functions <fi in 
the domain of H^ L : 

A, 

V n 4>(x) = —x(x)(j)(x) for x € 8Sl- 

For the formal definition of H^ L via quadratic forms and more background, 
see Mezincescu's original paper |25| as well as |231 124j . In particular, we 
immediately get the following important relations in the sense of the corre- 
sponding quadratic forms: 

H^ L < H Sl (24) 
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as well as 

n 

H^>@Hx lkiVk) , (25) 

k=l 

whenever the strip Sl is divided into disjoint strips Si k (yk) whose closures 
exhaust the closure of Sl- 

Proof of Proposition MIA Due to the form inequality (|24j) above it remains 
to prove the estimate for H^ L . 

Denoting the bottom eigenvalue of an operator H by E\ (H) (caution: here 
our notation differs from the one in |231 124j , where the second eigenvalue is 
denoted by E\(H)) we see that 

a(H^) n [E , E + L-T] + E 1 (H^) < E + L~\ 

Step 1. There exist b, K, f3 > such that 

P{E 1 {_H^ L ) <E + bL- 2 } < K ■ exp(-K ■ L dl ). (26) 

We use here the method from j^H]. Denote H(t,ui) := (Hq + tV u )^ L , and 
its first eigenvalue by E\(t,u). Since E\(t,uj) increases in t the event in 
(|26|) implies that E\(t,u) be small for all £ < 1 which in turn implies that 
E[(0, u) must be small. 

We infer from j24j . Theorem 3.25 that the gap between the first two 
eigenvalues satisfies 

E 2 (0,lo) -Ei(0,u) > const. L~ 2 . 
As in Lemma 2.3 this gives that 

\E 1 (t,uj)-{E + t-E[(0,uj))\<KL 2 -t 2 iorO<t <t-L- 2 . (27) 
Now assume that 

E X (H^) <E + bL~ 2 
for b > 0. From (j^J) we get that 

E[(0,u;) < c(b) 

with c(6) -► for 6 -► 0. 
On the other hand 

^(0,w) = (K^o,l|V<o,l) 

where ^o,z, is the normalized ground state of H ~,. Now, the boundary 
condition of H f, is defined so as to make sure that ipo is an eigenfunction; 
see the discussion in |24| . Therefore tpo,L = L _dl//2 V'o an d we get 

E'^u) = (K,Vo,l|V>o,l) 

= L-^Y,^) ■ [ U a (x)\Mx)\ 2 dx 

>L- dl il^-cu- ! |Vo(x)| 2 dx. 
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Since, by (B3), there are at least c±(L — rjj) dl elements of 1 n Sl~ Tu in 
R dl x A fix (0) and 



inf / \i/j (x) \ 2 dx > 

)£M d i xAb (0) JA rTr (xi.xi) 



(x 1 ,a!2)eM <i ixA J K_ L (0) yA rLr (xi,x 2 ) 

we arrive at 

^(0,w) >ci-|^-r £ »7«(w) (28) 

with c\ > and independent variables r/ Q running through an index set I±_ 
of cardinality at least c%L dl . If we now choose b > so small that < M 
where M is smaller than the mean of all the r/ a 's we get that: 

P{£l(iff) < £ + bL- 2 } < P{ Cl • -L £ ^(o;) < c (6)} 

<i4T ■exp(-^,|21|) 
< if • exp(-/?L dl ), 

by a standard large deviation estimate; see or |31| - Theorem 1.4. This 
finishes the proof of Step 1. 

Step 2. To deduce the desired bound from Step 1 we divide the strip Sl into 
disjoint strips Si k (yk) whose closures exhaust the closure of Sl and such that 

L~ 7 < b ■ ll 2 < 42 • L~ 7 , l k e 2N + 1 

which is possible for L large enough. 

Their number n is at most const, we know that 

so that 

rSr,\ ^ Z7_ l r H ^Tttf 17 t TJ Sl k^> k ' 



n 

< Y,y{Ei{Hx kiVk) ) < E + L-~<} 



k=l 
n 



<Y,HEi(Hx lk ' yk> )<E +b-if} 

k=l 

< n- K -exp(-pl k dl ) 

provided L is large enough. □ 

Remarks. (1) In cases where the operator H is ergodic, a stronger bound 
than (|23[) is provided in |241 Proposition 5.2]. Their bound is in terms of the 
integrated density of states for which |24| establishes Lifshits asymptotics. 
As we are only interested in localization properties here, the bound (|23|) 
suffices and allows to handle the non-ergodic random potentials defined in 
(B3) and (A3). 
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(2) We have established localization near the bottom of the spectrum 
for the random surface models considered in this section. If d\ = 1 one 
expects for physical reasons that the entire spectrum of H below inf a (Hp) 
(see (A4)) is localized. A corresponding result for lattice operators has been 
proven in ^Bj (in situations where Hp is the discrete Laplacian and di = 1). 
To show this for continuum models remains an open problem. 

5. Anderson models with displacement 

By considering the special case d\ = d, the results of the previous Section 
also cover "usual" Anderson models, sometimes also called alloy models. 
Note that in this case (B2) becomes trivial. Let us nevertheless state the 
assumptions and result again for this case, mainly because we want to point 
out below that the obtained bounds hold uniformly in the geometric param- 
eters describing the random potential. This will then be applied to models 
with random displacements. Here are the assumptions we rely upon: 

(Dl) Vq G Lf oc uni{ (M d ) is real-valued and periodic. 

(D2) The set X C Mr, where the random impurities are located, is uniformly 
discrete, i.e., inf{|a — f}\ : a ^ (3 G 1} =: rj > and uniformly dense, 
i.e., there exists Rj > such that A.r x (x) PiI ^ for every x G M. d . 

Theorem 11. Assume (Dl), (D2) and (A3). Then there exist 5 > 0, 
< s < 1, fi > and C < oo such that for I := [Eq,Eq + S], all open sets 
GcR d and x,yeR d 

sup n\\x x {H G -E-ie)- l Xy \\ s )<Ce-^~y\. (29) 

Eel,e>0 

In particular, the following consequences hold: 

(a) The spectrum of H G in I is almost surely pure point with exponentially 
decaying eigenf unctions. 

(b) There are fi\ > and C\ < oo such that for all x, y G Z d , 

E(sup||x.e- itHG P / (^ G )x y ||) < Cie-wl*-*!. (30) 

Here all the constants S,s,C, \i,G\, fj,i can be chosen to only depend on the 
potential through the parameters Vq, r/ max , M p , cjj, Cy, rjj, Ru, r%, Rj. 

To this end we first observe that (Dl), (D2) and (A3) imply (A4) with 
constants Ep, m and L* only depending on the listed parameters: 

Proposition 12. Assume (Dl), (D2) and (A3). Then there exist 

Ei = E x (Vq, r/ max , M p , cu, Cu, ru, Ru,r%, Rj) > E , 
m = m(V ,r] maiX , M p , cu, C v ,r v , Ru, r x , R%) G (0, 2) 

and L* = L*{. . .) such that 

(1) E F >E X . 

(2) For m d := 42 • d, all L > L* and x G Z d : 

P(<7(# Ai(x V)) n [£ ,£o + L~ m ] ± 0) < L" md . 
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Proof. First we show that (D2) implies that there exist cj, Cj and Lj de- 
pending only on n, Rj such that for all L > L%: 

c x -L d <#(lnA L (x))<C x -L d . (31) 

The upper bound follows from uniform discreteness: 

#(ln A L (x)) • \B rx/2 \ < \A L+rx/2 \ < (2L) d , 

provided L > rj/2. For the lower bound use uniform denseness: Divide 
Al(x) into disjoint boxes of side length Rj. If L > 2Rj there are at least 
(L/2Rj) d of them each of which contains at least one point from I. 

Now we can use the analysis of the preceding Section. Since the relevant 
quantities depend only on the indicated parameters, the assertions follow. 

□ 

With this uniform version of (A4) and the proofs provided in Sections |2 
and 01 we also get corresponding uniform versions of Theorems ^ and i.e. 
Theorem 1111 

As a specific application of the previous observation, we can start from 
an Anderson model as above and additionally vary the set X in a random 
way, as long as rj and Rj obey uniform upper and lower bounds. Instead 
of formulating the most general result in this direction we look at models 
that were introduced in ^0] and further studied in [3*2] . 

(D3) Let rfj, j € Z rf be independent random couplings, defined on a proba- 
bility space f2 with distribution pj and Uj as in (A3). 
(D4) Let Xj, j € Z d be independent random vectors of length at most | in 
denote the corresponding probability space by O,. 

Define 

H{u,u) := -A + Vq+Y, Vj(^)Uj{ - - j -Xj(^)). 

Corollary 13. Assume (Dl), (D3), (D4). Then, for H(u,u) as above there 
exist S > 0, < s < 1, fjL > and C < oo such that for I := [Eq,Eq + 6], all 
open sets G C M. d and x, y G M. d 

sup ^{\\x x {H G -E-\E)- x Xy \\ s )<Ce-^-y\. (32) 
Ee/,£>o 

In particular, the following consequences hold: 

(a) The spectrum of H G in I is almost surely pure point with exponentially 
decaying eigenf unctions. 

(b) There are \i > and C < oo such that for all x, y £ 1i d , 

EE(sup|| Xa; e- itHG P / (F G ) Xj/ ||) < Ce-^l. (33) 

teiR 

Proof. The corresponding inequality holds uniformly in u by what we proved 
above. □ 

Note that in this last Corollary we have not assumed that the random 
perturbations cover the whole space. In that respect our result provides 
substantial progress as compared to ^3 . 
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Appendix A. Some technical tools 

Here we collect some technical background which was used in Section 2 
above. All of this is known. We either provide references or, for convenience, 
in some cases sketch the proof. 

A.l. Combes-Thomas bounds. Proofs of the following improved Combes- 
Thomas bound can be found in (where it was first observed) and |30j . 
We state it here under assumptions which are sufficient for our applications. 
In particular, we assume d < 3, while the result holds in arbitrary dimension 
for a suitably modified class of potentials. As above, for an open G C M. d we 
denote by H the restriction of —A + V to L 2 (G) with Dirichlet boundary 
conditions. 

Proposition 14. Let d < 3, V G L 2 ocnni{ (R d ) with sup x \\VxAi( x )h < M - 
Let M > 1 and R > 0. Then there exist c\ = ci(M,R) and C2 = C2(M,R) 
such that the following conditions 

(i) Gcl d open, A, B C G, dist(A, B) =: 5 > 0, 

(ii) (r, s) C p(H G ) n (-R, R), E G (r, s) and r/ : = dist(£, (r, s) c ) > 0, 
imply the estimate 

sup \\xa{H g -E- ie)" 1 ^!! < ?l e -c2^v 1/2 s. (34) 

Note that the results in and [30] are stated for e = 0, but the proofs 
are easily adjusted to show that the bounds are uniform in the additional 
imaginary part. 

A. 2. Combes-Thomas bounds in Hilbert-Schmidt norm. A conse- 
quence of (|3*1)) is that \\xx{H^ G > — E — ie)~ 1 Xy\\ decays exponentially in 
\x — y\. Due to the restriction to d < 3 this is also true in Hilbert-Schmidt 
norm: 

Proposition 15. Let d < 3, V G ^i 2 OCiUnif (K d ), H = — A + V in L 2 (R d ) and 
L C (— oo, inf <t(H)) a compact interval. Then there exist C < oo and fi > 
such that 

sup \\ Xx (H G -E- ie^Xyhs < Ce~^ x -y\ (35) 

E&I,e>0 
GcM. d open 

for all x, y G M d . 

Proof. Let us sketch the proof by combining several well known facts. To 
this end, let S p denote the p-th Schatten class, i.e. the set of all bounded op- 
erators A such that \\A\\ P := (tr \A\ p yl p < oo. As d < 3, by Theorem B.9.3 
of |2E1 we have 

UAH - Ey^X < C! <oo (36) 

for each p > 3 and E < inf cr(H). The proof provided in [2Hj shows that 
C\ can be chosen uniform in x G M. d and E G /. In the sense of quadratic 
forms it holds that H G > H for each open G C R d , i.e. - Ey/ 2 (H G - 
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Ey 1/2 \\ < 1 for all E < mfa(H), see e.g. Section VI.2 of Thus 

\\ Xx {H G - E)-y% < \\ Xx (H - E)-V% \\(H - Ef'\H G - E)~W\\ 

< d < oo. (37) 
The Holder property of Schatten classes implies 

\\ Xx {H G -E^XyW^^Cl (38) 
uniformly in x, y G M. d , E G / and G C M. d open. From the resolvent identity 
Xx {H G -E-ie)- l Xy = x^-Ey'xy+iexA^-E-iey'i^-Er'xy 
we easily see that 

\\ Xx (H G -E-ie)- l Xy \\ p/2 <C 2 <oc (39) 

holds uniformly also in the additional parameter e € M.. By Proposition 1141 
we also have C3 < 00 and Hi > such that 

\\ Xx (H G -E- ie)- l x y \\ < C 3 e-«l^l (40) 

uniform in G, E £ I and e G R. As we may choose p/2 £ (3/2,2), (|35|) 
follows from (|39|) and (|40)1 by interpolation, more precisely from the fact that 
IHIhs = II -lb and PHI = tr |A| 2 = tr(\A\P/ 2 \A\ 2 -P/ 2 ) < p|| 2 -P/ 2 p||^. 

□ 

A. 3. A fractional-moment bound. The next result and its proof are 
found in [2] , where it played a central role in the extension of the fractional- 
moment method to Anderson-type random Schrodinger operators in the 
continuum. 

Recall that an operator A is called dissipative if Im(Aip, (p) > for all 
ip G E)(A). It is called maximally dissipative if it has no proper dissipative 
extension. Below we also use the notation | • | for Lebesgue measure in M. 2 . 

Proposition 16. There exists a universal constant C < 00 such that for 
every separable Hilbert space TL, every maximally dissipative operator A in 
TC with strictly positive imaginary part (i.e. Im(Aip,tp) > 5||<^|| 2 for some 
5 > and all ip G D(A)), for arbitrary Hilbert- Schmidt operators Mi, M 2 
in TC, for arbitrary bounded non-negative operators U\, U2 in TL, and for all 
t > the following holds: 

\{{vi,V2) 6 [0,1] 2 : \\M 1 ul /2 (A-v 1 U 1 -v 2 U 2 )~ 1 U 2 /2 M 2 \\ HS > t}\ 

<C||M 1 ||hs||M 2 ||hs- (41) 

The weak-Li-type bound 1)41(1 yields a fractional moment bound: 

Corollary 17. Let s G (0,1). Then for the constant C and operators A, 
Mi, M 2 , U\, U 2 as in Provosition \lb\ 

[ d Vl f dv^hU^iA-ViUi-v^^U^M^ 
Jo Jo 

<-^— HMxIlfjsHMsll^. (42) 
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This follows with layer-cake integration, which gives for the l.h.s. of Q42JI 

/ dui / d V2 ||...|| s < / \{(v uV2 ) € [0,1] 2 : ||...|| >t^ s }\dt. 

The integrand is bounded by the minimum of 1 and a bound following from 
(|4*T|) . Splitting the integral accordingly leads to (|i2~)) . 

Remarks. (1) The use of the interval [0, 1] as support of v\,V2 in Proposi- 
tion^] and Corollary 1171 is not essential. Using shifting and scaling it can 
be replaced by an arbitrary compact interval K, with constants becoming 
-RT-dependent. 

(2) In our applications maximally dissipative operators arise in the form 
A = — (S — E — ie) for self-adjoint operators S, with e > providing a 
strictly positive imaginary part. 

(3) Note that, as seen from the argument in .2:, a bound like Q42[) also 
holds in the "diagonal" case, i.e. for dv \\MU l l 2 (A - vU)' 1 ^ 2 M\\ s m . 
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